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Abstract 

The explosion of core-collapse supernova depends on a sequence of events taking place in less than a second 
in a region of a few hundred kilometers at the center of a supergiant star, after the stellar core approaches the 
Chandrasekhar mass and collapses into a proto-neutron star, and before a shock wave is launched across the 
stellar envelope. Theoretical efforts to understand stellar death focus on the mechanism which transforms the 
collapse into an explosion. Progress in understanding this mechanism is reviewed with particular attention 
to its asymmetric character. We highlight a series of successful studies connecting observations of supernova 
remnants and pulsars properties to the theory of core-collapse using numerical simulations. The encouraging 
results from first principles models in axisymmetric simulations is tempered by new puzzles in 3D. The 
diversity of explosion paths and the dependence on the pre-collapse stellar structure is stressed, as well 
as the need to gain a better understanding of hydrodynamical and MHD instabilities such as SASI and 
neutrino-driven convection. The shallow water analogy of shock dynamics is presented as a comparative 
system where buoyancy effects are absent. This dynamical system can be studied numerically and also 
experimentally with a water fountain. The potential of this complementary research tool for supernova 
theory is analyzed. We also review its potential for public outreach in science museums. 
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1 INTRODUCTION 

The explosive death of massive stars is a key ingre¬ 
dient in stellar evolution, stellar population synthesis 
and the chemical enrichment of galaxies. It defines the 
birth conditions of neutron stars which, if associated 
in a coalescing binary system, may be responsible for 
short GRBs and r-process nucleosynthesis. The shock 
wave launched in the interstellar medium during a su¬ 
pernova explosion accelerates cosmic rays and may con¬ 
tribute to triggering star formation. Understanding the 
mechanism of supernovae explosions has become a pri¬ 
ority since the detailed observation of SN1987A neutri¬ 
nos and the identification of its massive progenitor. It 
is still a theoretical challenge despite hundreds of new 
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events observed every year in distant galaxies, and the 
progress of computational power. Even if robust explo¬ 
sions are not obtained from first principles yet, numer¬ 
ical simulations are able to explore physical ideas in 
full 3D models, filling the gap between theoretical con¬ 
cepts and observations. The multiplication of theoreti¬ 
cal results may seem difficult to interpret in view of the 
multiplicity of physical assumptions, which range from 
simple adiabatic approximations of an ideal gas to ad¬ 
vanced modeling of neutrino transport and nuclear in¬ 
teractions in general relativity. Substantial progress has 
been achieved over the last ten years with an emphasis 
on the asymmetric character of the explosion, which we 
propose to summarize in the present review. 

Some important observational constraints regarding 
the asymmetric character of the explosion are recalled 
in Sect. |4.2[ The general theoretical framework of 
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neutrino-driven explosions is summarized in Sect. [3] 
In order to underline the nature of recent theoretical 
progress we have shown in Sect. [4] that a large set of 
modern physical ideas were already known before 2009 
and later confirmed by more advanced calculations. We 
analyze in Sect. [5] the new ideas which have driven su¬ 
pernova theory in the most recent years beyond the 
production of better tuned models. In particular the 
diversity of explosion scenarios has drawn attention to 
the pre-collapse stellar structure and correspondingly 
increased the size of the parameter space of initial con¬ 
ditions to be explored. This complexity calls for a deeper 
understanding of the physical processes. An experimen¬ 
tal fountain based on a shallow water analogy has been 
proposed to gain insight into one of the instabilities re¬ 
sponsible for the asymmetric character of the explosion. 
The potential of this new tool is analyzed in Sect. [6j 
both for theoretical research and for public outreach. 
The study of neutrino and gravitational wave diagnos¬ 
tics are left to the recent reviews of Janka et al. (2012), 
Janka (2012), Kotake et al. (2012), Burrows (2013) and 
Kotake (2013). Particular attention is paid to the many 
new results discovered over the last two years since these 
reviews. 


2 OBSERVATIONAL EVIDENCE FOR 
ASYMMETRIC EXPLOSIONS 

The observed distribution of pulsar velocities (e.g. 
Hobbs et al. 2005) has been a puzzle for more than 40 
years because they are much faster than massive stars 
(Gunn & Ostriker 1970). Average velocities of several 
hundreds of km/s cannot be explained by a residual or¬ 
bital velocity gained through the disruption of a binary 
system during the explosion. A pulsar kick is a natu¬ 
ral outcome of an asymmetric explosion process with a 


significant l = 1 component as seen in Sect. 4.2.1 The 
angle between the direction of the kick and the direc¬ 
tion of the rotation axis of the pulsar contains a very 
interesting constraint on the geometry of the explosion. 
The direction of the rotation axis can be accurately de¬ 
termined when a pulsar wind nebula is observed, but 
this sample is still small and lacks unambiguous cases 
where the kick is strong enough to neglect a possible 
orbital contribution to the kick (Ng & Romani 2004, 
Wang et al. 2006). Polarimetry suggests a strong corre¬ 
lation but cannot fully disentangle kick-spin alignment 
from orthogonality (Noutsos et al. 2012). This correla¬ 
tion is smeared out by the Galactic potential for pulsars 
older than 10 Myr (Noutsos et al. 2013). 

The asymmetric character of the explosion is also sug¬ 
gested by the sudden increase of polarization of the light 
observed from a type II-P supernova explosion after ^90 
days, at the moment of transition from the photospheric 
phase to the nebular phase (Leonard et al. 2006). This 
increase coincides with the moment when the innermost 


regions become visible. 

On a longer timescale, the shape of the inner ejecta 
currently seen in SN1987A also suggests an asymmetric 
explosion geometry which is not aligned with the large 
scale structure of the circumstellar medium (Larsson et 
al. 2013). Even three hundred years after the explosion, 
its asymmetric character can leave an imprint on the 
chemical composition of the ejecta. The spatial distri¬ 
bution of 44 Ti observed by the NuSTAR space telescope 
in Cassiopea A shows a global asymmetry, which seems 
to favor the region opposed to the direction of the com¬ 
pact object (Grefenstette et al. 2014). This one-sided 
asymmetry seems to confirm the theoretical prediction 
of Wongwathanarat et al. (2013) where a strong one¬ 
sided shock expansion is favorable to nucleosynthesis 
while the closer region in the opposite direction attracts 
the neutron star gravitationally. 

Direct constraints on the first second of the explosion 
are expected from the future detection of neutrinos (e.g. 
Wurrn et al. 2012, Muller et al. 2014a) and gravitational 
waves (Ott 2009, Kotake 2013, Muller et al. 2013) from 
a Galactic supernova. The future detection of the diffuse 
neutrino background will globally constrain both the 
physics of the explosion and the supernova rate (Bea- 
com 2010). For an individual Galactic event, the time 
variability of the neutrino signal from the IceCube ex¬ 
periment could directly measure the modulation of the 
emission due to a SASI oscillation of the shock (Lund 
et al. 2010, 2012, Tamborra et al. 2013, 2014). 


3 THEORETICAL FRAMEWORK 
3.1 Neutrino-driven explosion scenario 


The most promising framework to understand the ma¬ 
jority of supernova explosions was set by Bethe & Wil¬ 
son (1985). The delayed explosion mechanism driven 
by neutrino energy was initially described as a spher¬ 
ically symmetric process. As the mass of the central 
core approaches the Chandrasekhar mass (~ 1.4M so i), 
the pressure of degenerate relativistic electrons becomes 
insufficient to resist gravity. Electron pressure is further 
decreased as they are captured by protons and produce 
neutrons and neutrinos. A new equilibrium is reached 
between gravity and nuclear matter essentially made 
of neutrons. The stellar core collapses from a radius 
of ~ 1500km to this new equilibrium of a few tens of 
kilometers in less than half a second (Fig. [I]). The grav¬ 
itational energy gained from this collapse seems large 
enough to be able to account for the observed kinetic 
energy ~ 10 51 erg of the supernova ejecta: 


GMl 

fins 


2 x 10 53 erg 


/ 30knr\ / M ns \ 


V Rns ) V 1.5M so l/ 


( 1 ) 


A major difficulty of supernova theory is to explain how 
this gravitational energy is transferred to the stellar en- 
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Figure 1 . The neutrino-driven delayed explosion mechanism re¬ 
lies on the absorption of neutrinos by the dense post-shock gas. 


veloppe to reverse the infalling motion into an explo¬ 
sion. The infall of supersonic matter onto the surface of 
this proto-neutron star produces a deceleration shock 
which stalls at a radius of ~ 150km. The bounce of the 
free-falling matter Vg/2 ~ GM/r onto the core is far 
from elastic, because a significant fraction of the kinetic 
energy is absorbed into the dissociation of iron nuclei. 
This can be seen by comparing the kinetic energy of a 
free-falling nucleon and its binding energy 8.8MeV in 
the iron atom, neglecting special relativistic corrections 
for the sake of simplicity: 

\m n Vg 236km M 

8.8MeV r 1.5M sol ' ^ 

Stellar matter passing through the shock is thus dissoci¬ 
ated into free nucleons. The subsonic advection of these 
nucleons towards the surface of the proto-neutron star 
takes place in an intense flux of neutrinos which diffuse 
out of the proto-neutron star and carry away most of 
the gravitational energy gained during the contraction. 
The detection of neutrinos from SN1987A were instru¬ 
mental to confirm the premises of supernova theory set 
by Colgate & White (1966). 

The post-shock region is made of two successive regions 
defined by the direction of the reaction 

p + e O n + v, (3) 

where the relativistic velocity of the electrons corre¬ 
sponds to a Lorentz factor y e > (m n — m p )/m e ~ 2.5. 


Although most neutrinos emerge from the neutri- 
nosphere of the cooling neutron star, a fraction comes 
from the neutronization of accreted matter below the 
gain radius. This process decreases the entropy of neu¬ 
tron rich matter settling in a stably stratified manner. 
Between the shock surface and the gain radius, the dom¬ 
inance of neutrino absorption over neutrino emission in¬ 
creases the entropy of the matter. This gain region is 
crucial to the success of the explosion, which relies on 
the absorption of sufficient neutrino energy to revive the 
stalled shock. 

Other sources of energy have been considered to pro¬ 
duce an explosion, such as the rotational and the mag¬ 
netic energies. The rotation period of pulsars at birth 
defines a reference rotational energy E Iot deduced from 
the conservation of angular momentum 


-Grot 


2.8 x 10 5O erg 


( Mns \ 

\ 1.4M sol ) 


f Rn 


V 10km 
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10ms 



Spin periods P ns ~ 10 — 20ms at birth are considered 
by Heger et al. (2005) as plausible extrapolations from 
the observations of young pulsars, without excluding 
possible effects of binary interaction during the lifetime 
of the massive star (Sana et al. 2012, De Mink et al. 
2013, 2014) or a significant redistribution of angular 
momentum during or immediately after the pulsar 
birth. The relative inefficiency of known spin-down 
mechanisms (Ott et al. 2006) rules out rotationally- 
driven supernovae as the generic case. According 
to Eq. Q, the rotational energy could be a major 
contributor in the particular cases where the rotation 
rate of the stellar core is large enough to produce 
millisecond pulsars. 

Magnetic energy could also play an important role if 
the magnetic field within the core were strong enough. 
The main source of magnetic field amplification is the 
differential rotation within the core (e.g. Akiyama et 
al. 2003), which is only a fraction of the total rotational 
energy. 

Assuming that the rotational energy is too weak to 
be the dominant contributor of the most common 
explosions (i.e. P ns > 10ms), most studies of core 
collapse have focused on the challenge of producing 
an explosion without any rotation at all. We shall 
see in Sect. |4.3| that rotation could be an important 
ingredient of the explosion mechanism even if its total 
energy is modest. 

We can also expect magnetic effects to play a role 
in shaping the geometry of the explosion even if the 
magnetic energy is modest compared to the final 
kinetic energy of the ejecta (e.g. Guilet et al. 2011). 


Detailed numerical modeling solving the Boltzmann 
equation to describe neutrino transport and taking into 
account special and general relativistic effects reached 
the conclusion that the scenario proposed by Bethe & 
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Wilson (1985) is unable to power the spherically sym¬ 
metric explosion of massive stars (Liebendoerfer et al. 
2001) except for the lightest ones. The envelope of stars 
in the range 8-10 M so \ is light enough to allow for a 
spherical explosion of ~ 10 50 erg in about 100ms (Ki- 
taura et al. 2006). In a sense, the absence of systematic 
explosions in spherical symmetry is consistent with the 
observational evidence summarized in Sect. 14.21 Theo¬ 
retical efforts over the past two decades have updated 
the delayed neutrino-driven explosion scenario by tak¬ 
ing into account its multidimensional nature. 


SASI mechanism 



3.2 Hydrodynamical sources of asymmetry 


3.2.1 Neutrino-driven buoyancy 

In addition to the prompt convection associated with 
the deceleration of the shock as it stalls, the continued 
heating of the gas by neutrino absorption maintains a 
radial entropy gradient oriented in the same direction as 
gravity from the shock to the gain radius (Herant et al. 
1992, Janka & Muller 1996). Some gravitational energy 
can be gained if the flow is able to interchange high 
and low entropy layers. The Brunt-Vaisala frequency 
wbv characterizes the timescale of the fastest motions 
fed by buoyancy forces. Defining the entropy S' in a 
dimensionless manner, the Brunt-Vaisala frequency is 
expressed as follows: 


S = 


1 i P / P o 
7 ~ 1 ° S (p/Po) 7 ’ 


to 


2 

BV 


7 — 1 

- NS ■ V$. 

7 


( 5 ) 

( 6 ) 


This interchange can feed turbulent motions and push 
the shock further out, increasing the size of the gain 
region and diminishing the energy losses due to disso¬ 
ciation. Numerical simulations performed in the early 
2000’s were disappointing though (Buras et al. 2003), 
both because the duration of the simulation was lim¬ 
ited to a few hundred milliseconds and as some equa¬ 
torial symmetry was assumed due to limited compu¬ 
tational resources. The equatorial symmetry precluded 
the growth of global 1 = 1 modes. Foglizzo et al. (2006) 
pointed out that the negative sign of the entropy gradi¬ 
ent is not a sufficient criterion for the convective in¬ 
stability in the gain region because the interchange 
has to be fast enough to take place before the ad- 
vected gas reaches the gain radius. A criterion for lin¬ 
ear stability compares the advection timescale to the 
buoyancy timescale estimated from the Brunt-Vaisala 
growth rate: 


X = 


/•shock 


I^bvI 


dr 

kl 


< 3 


( 7 ) 


From this criterion one can anticipate that the strength 
and consequences of neutrino-driven buoyancy may 


Figure 2. SASI mechanism based on the coupling between acous¬ 
tic waves (wavy arrows) and advected perturbations (circular ar¬ 
rows) between the stalled shock and the proto neutron star. 


vary from one progenitor to another, depending on the 
radius of the stalled shock. 


3.2.2 Instability of the stationary shock: SASI 
Another source of symmetry breaking was discovered 
by Blondin et al. (2003) who studied a simplified setup 
where neutrino heating was neglected in order to avoid 
any possible confusion with neutrino-driven convection. 
The Standing Accretion Shock Instability (SASI) corre¬ 
sponds to a global oscillatory motion of the shock sur¬ 
face with a period comparable to the advection time 
from the shock to the neutron star surface. This global 
(l = 1) and oscillatory character are distinct features 
of the linear regime when perturbations grow exponen¬ 
tially with time. In contrast, the convective instability 
is non-oscillatory and dominated by smaller azimuthal 
scales (typically 1 = 5 — 6) comparable to the radial size 
of the gain region. 

The mechanism responsible for SASI relies on the in¬ 
teraction of acoustic perturbations and advected ones, 
such as entropy and vorticity perturbations (Galletti & 
Foglizzo 2005, Ohnishi 2006, Foglizzo et al. 2007, Scheck 
et al. 2008, Guilet & Foglizzo 2012). As illustrated in 
Fig-i a perturbation of the shock surface produces en¬ 
tropy and vorticity perturbations which are advected 
towards the proto-neutron star surface. Their deceler¬ 
ation in the hot region close to the neutron star pro¬ 
duces an acoustic feedback which propagates towards 
the shock surface, pushes it and regenerates new ad¬ 
vected perturbations with a larger amplitude than the 
first ones. This advective-acoustic cycle has been de¬ 
scribed analytically in a Cartesian geometry (Foglizzo 
2009, Sato et al. 2009). The same coupling mechanism 
is responsible for the instability of the bow shock in 
Bondi-Hoyle-Lyttleton accretion (Foglizzo 2001, 2002, 
Foglizzo et al. 2005). 
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4 WHAT DID WE KNOW BACK IN 2009? 

4.1 First successful explosions in numerical 
simulations from first principles 

The discovery of the global l = 1 oscillations of SASI 
emphasized the importance of performing numerical 
simulations over the full sphere. Using a ray by ray ap¬ 
proximation of neutrino transport and relativistic cor¬ 
rections to the gravitational potential, Marek & Janka 
(2009) obtained the first successful explosions of 11.2 
and 15 M so \ progenitor models from first principles. The 
explosion was more difficult to obtain with the more 
massive progenitor due to the intense ram pressure of 
infalling matter. It took 800ms of post-bounce evolu¬ 
tion (and three years of computation) for the axisym- 
metric shock to reach an explosion after a long phase 
of SASI oscillations. The easier explosion of 2D simu¬ 
lations compared to radial ones was analyzed by Mur¬ 
phy & Burrows (2008) who showed that post-shock gas 
trapped in convective patterns is exposed for a longer 
time to the neutrino flux, thus accumulating enough en¬ 
ergy and entropy in the gain region to push the shock 
outwards. The turbulent motions induced by hydrody- 
namical instabilities also contribute to the shock revival 
by building up a turbulent pressure in the post-shock 
region as observed by Burrows et al. (1995) and more 
recently stressed by Murphy et al. (2013) and Couch & 
Ott (2014). 


4.2 Instabilities in the supernova core as an 
explanation for observed asymmetries? 

4-2.1 The mystery of pulsar kicks potentially 
explained 

The Garching team was first to recognize the potential 
of global l = 1 deformations to explain pulsar velocities 
of several hundreds of kilometers per second (Scheck et 
al. 2004). The global conservation of linear momentum 
implies that the final momentum of the neutron star is 
equal in magnitude and opposed in direction to the total 
momentum of the ejecta. The asymmetric distribution 
of ejected matter has been studied in a series of axisym- 
metric simulations where the neutrino luminosity was 
adjusted to trigger an explosion. They showed that the 
gravitational interaction pulls the neutron star in the 
direction of the closest dense regions of post-shock gas, 
over a timescale of several seconds which can exceed the 
timescale of direct advection of momentum from the ac¬ 
creted gas. Their statistical study produced a shape of 
the distribution of pulsar velocities in agreement with 
observational constraints (Scheck et al. 2006). The mag¬ 
nitude of the kick is a stochastic variable, while the 
explosion energies and timescales were found to be de¬ 
terministic in this study. The magnitude of the kicks 
obtained has been subsequently confirmed with a dif¬ 


ferent axisymmetric code by Nordhaus et al. (2010a, 
2012) which directly followed the motion of the neu¬ 
tron star. This kick mechanism has been confirmed in 
3D simulations by Wongwathanarat et al. (2010, 2013) 
who used an axis-free Yin-Yang grid rather than spher¬ 
ical coordinates. 

4-2.2 Enhanced mixing triggered by asymmetric 
shock in SN1987A 

The light curve of SN1987A, together with the early 
emergence of X-ray and gamma rays suggested that 
significant mixing disrupted the onion like structure 
of the star (McCray 1993, Utrobin 2004, Fassia & 
Meikle 1999), probably triggered by the interaction of 
the shock with composition interfaces. The efficiency 
of this mixing is enhanced if the shape of the shock 
is non-spherical, as shown by Kifonidis et al. (2006). 
This study considered shock deformations dominated 
by the modes l = 1,2 as naturally produced by SASI 
during the first second after shock bounce. The sub¬ 
sequent hydrodynamical evolution of the shock across 
the stellar envelope up to the surface showed final iron 
group velocities up to 3300km/s, strong mixing at the 
H/He interface, and hydrogen mixed down to velocities 
of 500km/s. These features are closer to observations 
than those obtained by the same team from a spherical 
shock dominated by the smaller scale convective insta¬ 
bility (Kifonidis et al. 2003). The viability of this pro¬ 
cess in 3D was later confirmed by Hammer et al. (2010), 
who followed the evolution of the shock calculated in 3D 
by Scheck (2007), across the stellar enveloppe. They 
found that the development of mixing instabilities is 
more efficient in 3D and allows for higher clump ve¬ 
locities due to the different action of drag forces. The 
effect of the progenitor structure was recently explored 
by Wongwathanarat et al. (2014) with a series of red 
and blue supergiants of 15M so i and 20M so i. 

4.3 First hints of a dominant spiral mode in 
3D and the spin rate of neutron stars 

As early as 2007, the first 3D simulations of the SASI 
instability were performed in the simplest adiabatic ap¬ 
proximation (Blondin & Mezzacappa 2007). Instead of 
the sloshing oscillation 1 = 1, m = 0 ubiquitous in 2D 
axisymmetric simulations, the 3D simulations revealed 
a dominant spiral mode l = 1, m = ±1 even though no 
angular momentum was contained in the free falling 
stellar matter. This opened the possibility that a non¬ 
rotating stellar core could give birth to a rotating neu¬ 
tron star, with the opposite angular momentum carried 
away by the ejecta of the supernova. Blondin & Mezza¬ 
cappa (2007) estimated that the final spin period of the 
neutron star could be as short as 50ms. They checked 
that these dynamics were not an artifact of the inflow¬ 
ing inner boundary condition by performing equatorial 
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simulations with a hard surface and a cooling function 
(Blondin & Shaw 2007). 

Another spectacular result was obtained by Blondin & 
Mezzacappa (2007) by considering for the first time the 
3D development of SASI in a rotating progenitor. They 
observed that the spiral mode is preferentially triggered 
in the same direction as the progenitor rotation, and as 
a consequence the angular momentum received by the 
proto-neutron star is opposed to that of the parent star. 
The phase of SASI instability may thus decelerate a 
neutron star to a lower spin than naively deduced from 
the conservation of angular momentum. They showed 
that a rotating massive star may even give birth to a 
counter rotating neutron star. These results were con¬ 
firmed by Yamasaki & Foglizzo (2008) using a pertur¬ 
bative approach in a cylindrical geometry. In particular 
the growth rate of SASI was shown to depend linearly 
on the angular momentum of the progenitor. Rotation is 
able to destabilize the prograde mode of SASI through 
the difference of rotation rates between the shock radius 
and the inner radius, even when the rotation rate is so 
small that the centrifugal force is negligible. 

The existence of a physical mechanism capable of mod¬ 
ifying the angular momentum of the neutron star dur¬ 
ing the collapse revived the debate about the compari¬ 
son between the angular momentum in the stellar core 
inferred from stellar evolution and the angular momen¬ 
tum of pulsars at birth inferred from observations. The 
study of Heger et al. (2005) had concluded that the spin 
rates of 10-15ms of the most common pulsars were com¬ 
patible with stellar evolution with plausible magnetic 
field strengths deduced from dynamo action, without 
the need to invoke any magnetic braking during or after 
the explosion. Alternate scenarios are possible though, 
since the transport of angular momentum in the 1-D 
codes of stellar evolution relies on uncertain prescrip¬ 
tions based on debated dynamo action (Spruit 2002, 
Zahn et al. 2007). In particular the distribution of an¬ 
gular momentum could also be affected by the action of 
internal gravity waves (e.g. Lee & Saio 1993, Talon & 
Charbonnel 2003, Pantillon et al. 2007, Lee et al. 2014). 

4.4 The uncertainties of a non-axisymmetric 
explosion scenario 

The non-axisymmetric shock geometry observed in the 
first 3D simulations by Blondin & Mezzacappa (2007) 
cast doubts on all the conclusions drawn in the axisym- 
metric hypothesis and urged new 3D simulations incor¬ 
porating non-adiabatic processes. Would the third di¬ 
mension be a key ingredient for robust successful explo¬ 
sions? Contrasting with the ample deformations of the 
shock found in the idealized simulations of Blondin & 
Mezzacappa (2007), the first non-adiabatic simulations 
by Iwakami et al. (2008) suggested that the asymme¬ 
tries induced by SASI may be weaker in 3D than in 2D. 


Back in 2009 it was very unclear whether the pulsar 
kick would reach a comparable magnitude in 3D and in 
2D, and whether the efficiency of mixing induced by the 
shock propagation in the stellar enveloppe would com¬ 
pare as favorably to SN1987A in 3D simulations as in 
2D. As already noted in Sect. |4.2[ these two questions 
later received positive answers. Launching an explosion 
from first principles, however, turned out to be even 
more challenging in 3D than in 2D. 

5 WHAT IS NEW SINCE 2009? 

The progress of supernova theory in the last 5 years 
has gone much further than confirming previous results 
with improved computational power. New constraints 
have been obtained from the determination of neutron 
star masses. New ideas have emerged concerning the 
diversity of explosion paths and their sensitivity to the 
structure of the progenitor above the iron core includ¬ 
ing its angular momentum, magnetic field and initial 
asymmetries. 

5.1 A better constrained equation of state 

The equation of state of neutron rich matter at nuclear 
densities is an important ingredient of the dynamics of 
core-collapse, which sets the radius of the proto-neutron 
star and thus the depth of the gravitational potential. 
Most simulations have used a parametrized equation 
of state provided by Lattimer & Swesty (1991), Hillc- 
brandt & Wolff (1985) or Shen et al. (1998). The de¬ 
scription of matter properties at nuclear densities is 
extrapolated from the properties of nuclei available in 
terrestrial experiments, as well as measurements of the 
mass and to a lesser extent the radius of neutron stars. 
The measurement of the mass of a neutron star with 
a mass of 1.97 ± 0.04M so i by Demorest et al. (2010) is 
accurate enough to rule out a series of alternate mod¬ 
els and significantly reduces the uncertainty associated 
to the equation of state. A second neutron star in this 
mass range has been discovered by Antoniadis et al. 
(2013), with a mass of 2.01 ± 0.04M so i for the pulsar 
PSR J0348+0432. Exotic forms involving a transition 
to quark matter (Sagert et al. 2009, Fischer et al. 2011) 
can produce significant dynamical consequences dur¬ 
ing core collapse, but involve phenomenological mod¬ 
els with parameters adjusted to reproduce the observed 
massive neutron stars. In particular, the properties of 
quark matter at the relevant densities are neither acces¬ 
sible to experiments nor to any theoretical development 
from first principles. 

Staying in the classical framework, some improvements 
are necessary to account for the latest advances in nu¬ 
clear physics. Several statistical models have been de¬ 
veloped accounting for the entire distribution of nuclei 
(Hempel & Schaffner-Bielich 2010, Raduta & Gulminelli 
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2010, Buyukcizmeci et al. 2013) in contrast to the ap¬ 
proximation of a single representative heavy nucleus 
and a-particles within the classical EoS (Lattimer & 
Swesty 1991, Shen et al. 1998). Except in some small 
density and temperature regions, the effects on global 
thermodynamics remains small and the effect on core 
collapse simulations is modest (Steiner et al. 2013). The 
effect on the neutrino interactions and spectra could be 
more important (Arcones et al. 2008). 

For progenitors with masses above ~ 25M so i, the tem¬ 
peratures and densities reached during core collapse can 
become so high that a traditional description in terms 
of electrons, nuclei, and nucleons is no longer adequate. 
Several EoS have been developed incorporating addi¬ 
tional particles, such as hyperons and pions (Oertel 
et al. 2012, Gulminelli et al. 2013, Banik et al. 2014). 
Among the consequences, the time necessary to collapse 
to a black hole is reduced (Peres et al. 2013). 

5.2 Confirmed axisymmetric explosions from 
first principles for a wider set of 
progenitors 

A series of paper by Suwa et al. (2010), Muller et al. 
(2012a, 2012b, 2013), and Bruenn et al. (2013) con¬ 
firmed the viability of the neutrino-driven explosion 
mechanism in the axisymmetric hypothesis. Successful 
explosions have been obtained from first principles for 
progenitor masses ranging from 8.1 to 27M so i, both by 
the MPA group and the Oak Ridge group. One may 
regret that the explosion times and energies measured 
by these two groups are not yet fully compatible with 
each other, and neutrino transport still relies on some 
numerical approximations, but these recent studies are 
still very encouraging steps towards a consensus. Both 
groups use a ray-by-ray method which solves the Boltz¬ 
mann equation in each radial direction assuming an ax¬ 
isymmetric distribution of neutrinos along each radial 
ray. ft is interesting to note that the Princeton group 
did not obtain an explosion using a different code where 
neutrino transport is approximated with a multi-group, 
flux-limited diffusion method (Dolence et al. 2014). 
Although these explosions are an immense success in 
view of the decades of failed numerical attempts, the 
mechanism cannot be qualified as robust enough yet, 
since the final explosion energy seems significantly too 
low compared to observations. It is not clear whether 
the solution to this puzzle will come from improving the 
modeling of progenitor models, or including 3D effects 
from rotation and magnetic fields, or something else. 

5.3 The diversity of explosion paths 

Despite the universality of the Chandrasekhar mass 
defining the initial conditions of core collapse, the past 
five years have taught us that the explosion process is 


sensitive to the radial structure of the envelope sur¬ 
rounding the central core. 

5.3.1 The parameter y and the relative roles of 
neutrino-driven convection vs SASI 

The numerical simulations performed at MPA and 
ORNL confirmed the diversity of hydrodynamical pro¬ 
cesses ruling the multidimensional evolution during 
the last hundreds of milliseconds before the explosion. 
Neutrino-driven buoyancy seems to govern the post¬ 
shock dynamics of the progenitor of 11.2 M so \. The 
global oscillations of SASI are clearly dominant for the 
27 M so i progenitor (Muller et al. 2012b), while both 
instabilities seem to be entangled in the evolution the 
15M so i progenitor. Such differences can be understood 
from the comparison of the advection time and the 
buoyancy times in these systems. A short advection 
time is both favorable to SASI (Foglizzo et al. 2007, 
Scheck et al. 2008), and stabilizing for neutrino-driven 
convection (Foglizzo et al. 2006). By adjusting the neu¬ 
trino luminosity and the rate of dissociation at the 
shock, Fernandez et al. (2014) were able to characterize 
each instability separately and compare the properties 
of the buoyant bubbles and turbulence induced before 
the explosion. 

Besides a few examples in the parameter space, we do 
not have a clear picture of the dependence of the y pa¬ 
rameter on the main sequence mass yet. It should be 
noted that the development of SASI is not a sufficient 
criterion for a successful explosion. The axisymmetric 
simulation of a 25M so i progenitor is a clear example of 
strong SASI activity without an explosion (Hanke et al. 
2013). 

5.3.2 The compactness parameter £ 2.5 and black 
hole vs neutron star formation 

A systematic study performed by Ugliano et al. (2012) 
used some educated prescription to account for multi¬ 
dimensional evolution using ID simulations over a wide 
range of progenitor masses from 10 to 40 solar masses. 
The outcome of the core collapse can either be a super¬ 
nova explosion with the formation of a neutron star, a 
failed supernova with the formation of a black hole, or 
a supernova with the delayed formation of a black hole 
through the fall back of enough ejecta. The diversity 
of mass loss processes during the life of a star results 
in a non-monotonous relation between its mass at the 
moment of core-collapse and its mass on the main se¬ 
quence. In addition to stressing this labeling confusion, 
the study of Ugliano et al. (2012) showed that the pro¬ 
duction of a neutron star or a black is also very variable 
for progenitors more massive than 15 M so \. An interest¬ 
ing indicator is based on the compactness parameter 
£ 2.5 defined by O’Connor & Ott (2011), measuring the 
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radius of the innermost 2.5M so i: 

t _ M/M S oi 

?2 ' 5 ~ i?(M)/1000km' [> 

Although this indicator is not strictly deterministic, a 
threshold £ 2 .5 > 0-3 seems very favorable to black hole 
formation in the calculations of Ugliano et al. (2012). 
From a series of 101 axisymmetric simulations from 10.8 
to 75M so i with IDS A neutrino transport, Nakamura 
et al. (2014a) found that a high compactness param¬ 
eter corresponds to a later explosion time, a stronger 
neutrino luminosity and explosion energy, and a higher 
Nickel yield. 

5.3.3 The possible influence of pre-collapse 
convective asymmetries 

The structure of the progenitor may also influence the 
asymmetry of the explosion mechanism through the 
convective inhomogeneities associated with thermonu¬ 
clear burning in the Silicon and Oxygen shells (Arnett & 
Meakin 2011). The quantitative impact of this process 
has been estimated by Couch & Ott (2013), who showed 
an example where these asymmetries had an effect com¬ 
parable to a 2% increase of the neutrino luminosity and 
could be enough to turn a failed explosion into a success¬ 
ful one. Elaborating in this direction, the recent study 
of Muller & Janka (2014b) pointed out that the most 
efficient effect of pre-collapse asymmetries comes from 
inhomogeneities with the largest angular scale l = 1,2. 

5.4 The unexpected difficulties of 3D 
explosions 

5-4-1 3D explosions are not easier than 2D 
Using a simplified model of neutrino interactions with 
cooling and heating functions, Nordhaus et al. (2010b) 
had claimed that the neutrino luminosity needed to ob¬ 
tain an explosion was lower in 3D than in 2D, in a 
comparable proportion as measured by Murphy & Bur¬ 
rows (2008) between 2D and ID. This raised the hope 
that the weak explosions obtained in 2D axisymmet¬ 
ric simulations could become more robust once com¬ 
puter ressources would allow first principle calculations 
including neutrino transport in 3D. A careful check 
by Hanke et al. (2012) and Couch (2013) showed that 
the neutrino luminosity threshold in the idealized setup 
studied by Nordhaus et al. (2010b) is actually very simi¬ 
lar in 3D and 2D. The misleading results of Nordhaus et 
al. (2010b) happened to be due to some numerical errors 
in the treatment of gravity in their CASTRO code, later 
corrected and acknowledged by Burrows et al. (2012). 
Despite these corrections, the convergence between the 
different codes is not yet fully satisfactory (Dolence et 
al. 2013). 

Using the IDSA approximation of neutrino transport, 
Takiwaki et al. (2014) studied the evolution of a 11.2 


M so i progenitor and observed that the explosion time is 
shorter and the explosion is more vigorous in 2D than 
in 3D. Similar conclusions were reached by Couch & 
O’Connor (2014) with a multi species leakage scheme 
for progenitors of 15 and 27 M so \. Let us note that the 
stochastic nature of the dynamical evolution and the 
prohibitive cost of 3D simulations make it difficult to 
achieve numerical convergence. Beside the possible in¬ 
fluence of numerical artifacts (e.g. Ott et al. 2013 com¬ 
pared to Abdikamalov et al. 2014), part of the difference 
between 2D and 3D simulations could be related to a 
more efficient energy cascade from large scales to small 
scales in 3D than in 2D (Hanke et al. 2012, Couch 2013, 
Couch & O’Connor 2014). 


5-4-2 SASI does exist in 3D models of core collapse 
According to the mechanism proposed by Guilet et al. 
(2010) and checked on the axisymmetric simulations of 
Fernandez & Thompson (2009), the saturation ampli¬ 
tude of SASI oscillations is not expected to be very dif¬ 
ferent in 2D and 3D. The waves of entropy and vorticity 
created by the oscillations of the shock are expected to 
lose their large scale 1 = 1 coherence when disrupted 
by the parasitic growth of Rayleigh-Taylor and Kelvin- 
Helmholtz instabilities. This description, however, ig¬ 
nores the interaction of SASI with small scale turbu¬ 
lence, whose properties may differ in 2D and 3D. Even 
though the first 3D simulations performed by Blondin 
& Mezzacappa (2007) in an adiabatic approximation 
showed well defined spiral shock patterns on a large 
scale, subsequent 3D simulations with less idealized se¬ 
tups suggested weaker amplitudes than in 2D (Iwakami 
et al. 2008, Burrows et al. 2012). The clear dominance of 
non-oscillatory convective motions in some simulations 
even led Burrows et al. (2012) to draw general conclu¬ 
sions such that SASI would be an artefact of simplified 
physics which does not exist in realistic 2D simulations, 
and even less in 3D. This reasoning missed the diver¬ 
sity of explosion paths later explored in a parametric 
manner by Fernandez et al. (2014) in 2D or Iwakami et 
al. (2014a) in 3D, and exemplified by the axisymmetric 
simulation by Muller et al. (2012) of the 27 M so \ pro¬ 
genitor. The evolution of this progenitor was computed 
in 3D from first principles by Hanke et al. (2013) who 
confirmed the presence of the SASI mode with an ampli¬ 
tude which can exceed even the amplitude observed in 
2D. Nevertheless, the explosion did not take place dur¬ 
ing the first 400ms of evolution, whereas it did explode 
in 2D on this timescale. It is interesting to note the 
sensitivity of the dynamics to small differences in the 
equation of state and possibly GR effects between the 
axisymmetric results of Hanke et al. (2013) and those 
of Muller et al. (2012), where less than 200ms were suf¬ 
ficient for an explosion. 
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5.5 Rotation and magnetic fields 

The angular momentum of the stellar core and its mag¬ 
netic field are two physical ingredients which definitely 
affect the birth properties of the neutron star and may 
also influence the explosion mechanism. 

5.5.1 SASI effect on the pulsar spin 

The spectacular conclusions regarding the pulsar spin 
obtained by Blondin & Mezzacappa (2007) with a 3D 
adiabatic setup were naturally followed by more realis¬ 
tic 3D simulations. Fernandez (2010) was able to con¬ 
firm the angular momentum budget when the progeni¬ 
tor does not rotate, with a setup similar to Blondin & 
Shaw (2007) but in 3D. By contrast, the spiral mode 
did not appear easily in the non-rotating simulation 
of Iwakami et al. (2009) who considered a stationary 
accretion flow incorporating the buoyancy effects due 
to neutrino heating and a more realistic equation of 
state. Significantly slower pulsar spin periods of the or¬ 
der of 500-1000ms were estimated by Wongwathanarat 
et al. (2010) who considered the collapse of a 15M so i 
progenitor with neutrino cooling and heating, self grav¬ 
ity and general relativistic corrections of the monopole, 
and a Yin-Yang grid. Similar skeptical conclusions were 
reached by Rantsiou et al. (2011), who pointed out that 
the mass cut during the explosion process might not co¬ 
incide with the separation of positive and negative an¬ 
gular momentum produced by SASI. A larger sample 
of non-rotating progenitors was considered by Wong¬ 
wathanarat et al. (2013), including a 20M so i progenitor, 
and measured final spin periods of 100ms to 8000ms. 
They noted that the kick and spin do not show any obvi¬ 
ous correlation regarding their magnitude or direction. 
Their analysis clarified an important difference between 
the kick and spin mechanisms: the spin up takes place 
during the phase of accretion onto the neutron star sur¬ 
face because it is dominated by the direct accretion of 
angular momentum. By contrast, the kick can grow on 
a longer time scale, even when accretion has stopped, 
through the action of the gravitational force. Using an¬ 
alytic arguments, Guilet & Fernandez (2014) estimated 
the maximum spin of a neutron star born from a non¬ 
rotating progenitor by measuring the amount of angular 
momentum which can be stored in the saturated spiral 
mode of SASI. They obtained spin periods compatible 
to those measured in the simulations of Wongwatha¬ 
narat et al. (2013). 

5.5.2 Rotation effects on the shock dynamics 

The effect of a moderate rotation on the dynamics of 
the collapse has received little attention so far, despite 
the fact that interesting correlations between the kick 
and spin directions would be expected in this regime. 
Indeed, the growth of a prograde spiral SASI mode is ex¬ 
pected from the early numerical simulations of Blondin 


& Mezzacappa (2007), and the linear stability analysis 
of Yamasaki & Foglizzo (2008) indicated that this ef¬ 
fect of differential rotation can be significant even if the 
centrifugal force is weak. The 3D simulations of sta¬ 
tionary accretion by Iwakami et al. (2009) confirmed 
that the amplitude of the prograde spiral mode of SASI 
is enhanced by rotation. Iwakami et al. (2014b) grad¬ 
ually increased the angular momentum of the infalling 
gas to measure the influence of rotation on the explo¬ 
sion threshold. According to the few points in their 
Fig. 1, the change in the threshold of neutrino lumi¬ 
nosity required to obtain an explosion seems to be a 
linear function of the angular momentum with a lumi¬ 
nosity variation of ~ 10% for an angular momentum of 
the order of ~ 5 x 10 15 cm 2 s _1 . More simulations would 
be necessary to check whether this effect really scales 
linearly with the rotation rate as expected. According 
to a tentative linear extrapolation, rotation effects with 
an angular momentum of ~ 10 15 cm 2 s _1 could induce 
a 2% shift of the neutrino luminosity, i.e. a similar 
order-of-magnitude effects as the pre-collapse asymme¬ 
tries described by Couch & Ott (2013). For reference, 
10 15 cm 2 s -1 is the equatorial angular momentum at the 
surface of a pulsar with a radius of 10km and a spin pe¬ 
riod of 6ms. The uncertain effect of mixing instabilities 
inside the neutron star remains to be evaluated. Com¬ 
pared to the models of stellar evolution by Heger et al. 
(2005), a value of 10 15 cnr 2 s _1 is twice as large as the es¬ 
timated angular momentum at 2000km with magnetic 
fields (~ 5 x 10 14 cm 2 s _1 ), and 20 times smaller than 
without magnetic fields (~2x 10 16 cm 2 s _1 ). 
Nakamura et al. (2014b) considered the collapse of the 
stellar core of a 15M so i progenitor with a shell type 
rotation profile: f2(r) = fl 0 /(l + r 2 /R 2 ) with R 0 = 2 x 
10 8 cm and Qq = 0, 0.17T, and 0.57T rad s _1 correspond¬ 
ing to an initial angular momentum at 2000km of 0, 
6 x 10 15 cm 2 s _1 , and 3 x 10 16 cm 2 s _1 . They measured a 
reduction of the threshold of neutrino luminosity of the 
order of 10% for Uo = 0.1-7T, which is consistent with the 
results of Iwakami et al. (2014b). They emphasized that 
the direction of the explosion is preferentially perpen¬ 
dicular to the spin axis but did not take into account 
the asymmetry of neutrino emission induced by rota¬ 
tion, which may favour the axial direction. 

5.5.3 The growth of magnetic energy without 
differential rotation 

Endeve et al. (2010, 2012) used an adiabatic setup to 
show that even when the initial rotation of the core is 
neglected, the magnetic energy can grow from the tur¬ 
bulent motions induced by SASI to reach 10 14 G at the 
surface of the neutron star. The growth of this magnetic 
field, however, did not significantly affect the dynam¬ 
ics of the shock. Obergaulinger et al. (2014) addressed 
the same question in a more realistic setup including 
neutrino-driven convection and a Ml approximation of 
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neutrino transport. They confirmed the lack of signif¬ 
icant effects on the shock evolution except if the field 
is initially very strong. They found evidence for an ac¬ 
cumulation of Alfven waves at the Alfven point as de¬ 
scribed by Guilet et al. (2011), but did not find a signif¬ 
icant contribution to the field amplification or heating. 

5.5.4 Magnetic field amplification with fast rotation 
Fast differential rotation provides a large energy reser¬ 
voir for magnetic field amplification, such that magnetic 
effects cannot be neglected beyond a certain angular 
momentum. The shearing of a poloidal magnetic held 
into a toroidal one provides a linear amplification, which 
would be relevant only if the initial poloidal magnetic 
held were large enough. If the initial poloidal magnetic 
held is weak, the exponential growth of the magnetoro- 
tational instability (MRI) is a more promising mech¬ 
anism of magnetic held amplihcation (Akiyama et al. 
2003). A meaningful description of this process requires 
3D simulations with a very high resolution, due to the 
very short wavelength of the MRI growing on a weak ini¬ 
tial held. The prohibitive computing time required has 
led researchers to use axisymmetric simulations with 
artificial assumptions such as a strong initial poloidal 
held mimicking the outcome of some amplihcation pro¬ 
cesses (e.g. Moiseenko et al. 2006, Burrows et al. 2007, 
Takiwaki et al. 2009, Takiwaki & Kotake 2011). In this 
framework, they obtained powerful magnetorotational 
explosions, with jets launched along the polar axis due 
to the winding of the initial poloidal magnetic into a 
very strong (> 10 15 G) toroidal held. The progress of 
computational power has opened new perspectives to 
address this problem in 3D. Winteler et al. (2012) sug¬ 
gested that such magnetorotational explosions could be 
an important source of r-process elements. The 3D sim¬ 
ulations of Mosta et al. (2014) highlighted the role of a 
non-axisymmetric instability of the toroidal magnetic 
held which can break the jet structure and prevent 
such magnetorotational explosions, a non-runaway ex¬ 
pansion of the shock being observed instead. 

In parallel progress is being made in understanding the 
magnetic held amplihcation due to the MRI in the pro¬ 
toneutron star. Buoyancy driven by radial gradients of 
entropy and lepton fraction has a strong impact on the 
dynamics, though the thermal and lepton number dif¬ 
fusion due to neutrinos allows for fast MRI growth, as 
shown by the linear analysis of Masada et al. (2006, 
2007). The study of the non-linear phase of the MRI 
requires 3D simulations, which have been undertaken 
in local and semi-global models that describe a small 
fraction of the protoneutron star. Taking into account 
global gradients responsible for buoyancy has, however, 
proven difficult in this framework because of boundary 
effects (Obergaulinger et al. 2009) or coarse radial reso¬ 
lution (Masada et al. 2014). A local model in the Boussi- 
nesq approximation allowed Guilet & Muller (2015) to 



Figure 3. The shape of the rotating hydraulic jump observed in 
the non-linear regime in the SWASI experiment (right) is sim¬ 
ilar to the shape observed in the shallow water approximation 
(left) and the shape of the SASI in the numerical simulations of 
cylindrical gas accretion with local neutrino cooling (top). 


estimate that a stable stratification decreases the effi¬ 
ciency of magnetic field amplification only slightly, but 
significantly impacts the structure of the magnetic field. 
Neutrino radiation, neglected in most numerical simu¬ 
lations, can slow down the MRI growth through an ef¬ 
fective neutrino viscosity deep inside the proto-neutron 
star or through a neutrino drag near its surface (Guilet 
et al. 2014). Further local simulations will be useful to 
study the MRI in this new growth regime, and ulti¬ 
mately global simulations (only achieved in 2D so far 
(Sawai et al. 2013, Sawai & Yamada 2014)) will be 
needed to assess the magnetic field geometry and its 
impact on the explosion. 

6 AN EXPERIMENTAL APPROACH TO 
SUPERNOVA DYNAMICS 

In order to make the SASI phenomenon more intuitive, 
Foglizzo et al. (2012) designed a shallow water analogue 
of the gas motion in the equatorial plane of the stel¬ 
lar core. Surface gravity waves in water play the role 
of acoustic waves in the stellar gas and participate in 
an unstable cycle with vorticity perturbations. Solving 
numerically the 2D system of shallow water equations 
revealed the similarity between SASI and its shallow 
water analogue. The absence of buoyancy effects in this 
shallow water formulation is theoretically instructive, 
as will be shown further. The first laboratory analogue 
of SASI dynamics was obtained using a water fountain 
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with a diameter of 60cm. This new research tool can 
help us become more familiar with SASI, in the same 
manner as we are familiar with the flapping of a flag or 
the convection of a fluid heated from below. 


6.1 From astrophysical complexity to the 
simplicity of a water fountain 

6.1.1 From gas to water dynamics: the shallow 
water analogy 

In the simple fountain built by Foglizzo et al. (2012), the 
dynamics of water looks very similar to the gas dynam¬ 
ics induced by SASI in the equatorial plane of the star, 
on a scale one million times smaller and an oscillation 
period one hundred times slower than in the astrophys¬ 
ical flow. On the one hand, the fury of extreme gravity 
accelerating the stellar gas to one tenth of the speed of 
light in less than half a second. Iron nuclei broken apart 
upon compression by a shock wave which moves back 
and forth with a 30ms period. On the other hand, wa¬ 
ter flowing at room pressure and temperature in a meter 
size fountain, giving rise to the back and forth motions 
of a surface wave with a 3 second period. The striking 
resemblance illustrated in Fig. [3]is rooted in the similar¬ 
ity between the Euler equations describing an inviscid 
adiabatic gas and the Saint Venant equations describ¬ 
ing shallow water. The conservation of mass, momen¬ 
tum and energy for a gas with an adiabatic index 7 , 
density p , pressure P , entropy S, velocity v, and sound 
speed c 2 = qP/p is written as follows: 


dv 

—+ (Vx»)x« + V 


^ + V • (pv) = 0, (9) 


- 1 -— 

2 7 -I 

dS 


+ $ 


= <£v#L0) 


— +V .S7S = 0. (11) 


The pressure gradient VP in the Euler equation 
( 10 ) has been decomposed into a enthalpy gradient 


pV(Cg /(7 — 1 )) and an entropy gradient pc^VS 1 . 

The Saint Venant set of equations describing the shal¬ 
low layer of water with depth H flowing with a velocity 
v over a surface z = H${r) is the following: 


dH 

ht 


+ V • (Hv) = 0, 


dv 

dt 


+ (V x v) x 'C + V 



= 0 . 


( 12 ) 

(13) 


The local gravity g = 981cm s -2 in the laboratory in- 
fuences the speed c w of surface waves in the experiment 
and the effective potential $ = gH $, which mimicks the 
Newtonian gravity of the central object with a hyper- 
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bolic shape: 


Cw = gH, 

H$(r) = — (5 - 6cm ) 2 


(14) 

(15) 


The Froude number defined by Fr = |u|/Cw plays the 
same role as the Mach number M. = \v\/c s in a gas. 
The pressure variation A P induced by the weight of 
water A P = p w g(z — P$(r)) has a mean value noted 
P = p w gH/2 related to the local speed of surface waves 
c w - 


p = 



(16) 


This relation is analogous to the relation P = pc 2 /q for 
a gas with an adiabatic index 7 = 2. 


6.1.2 The absence of buoyancy effects in shallow 
water 

Besides the particular value of the adiabatic index, the 
main difference between the two sets of equations comes 
from the absence of buoyancy effects in shallow wa¬ 
ter. This limitation precludes the study of buoyancy 
driven instabilities such as neutrino driven convection 
and their interesting interaction with SASI. This lim¬ 
itation also ignores the formation of buoyant entropy 
bubbles from shock oscillations driven by SASI (e.g. Fer¬ 
nandez et al. 2014), and the dominant role of buoyancy 
driven parasitic instabilities to saturate the amplitude 
of SASI (Guilet et al. 2010). Even in the inviscid limit 
the dynamical system described by the shallow water 
equations can only describe the interaction of surface 
waves and vorticity perturbations, whereas the Euler 
equations describe the interaction of acoustic waves, en¬ 
tropy and vorticity perturbations. In a positive sense, 
the shallow water dynamical system offers an opportu¬ 
nity to assess the importance of buoyancy effects on the 
dynamics of the shock. 


6.1.3 Inner boundary condition 

The inner boundary condition in the experiment faces 
the same difficulty as the adiabatic simulation of 
Blondin & Mezzacappa (2007): without any non- 
adiabatic process to shrink the volume of the accreted 
fluid on the surface of the neutron star, the fluid has to 
be continuously extracted through the inner boundary. 
This is done in the experiment by allowing water to spill 
over the upper edge of the inner cylinder (Fig. [4]). As 
long as the radius R ns of the cylinder is large enough to 
allow the free fall of water inside it, this inner bound¬ 
ary condition can be translated into the condition that 
the radial component of the flow velocity |u ns | reaches 
the local wave velocity q, (P ns ) associated to the depth 
Hedge of water above the upper edge of the inner cylin- 
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Figure 4. Water is injected inward from a circular slit, in a uni¬ 
form and stationary manner. A circular hydraulic jump is pro¬ 
duced by the vertical surface of the cylinder at the center. Water 
is evacuated by spilling over the upper edge of this cylinder. 


der. Denoting by Q the flow rate, 

2TrR ns H edge \v ns \ = Q, (17) 

^ns = Sledge- ( 18 ) 

6.1.4 SWASI, a shallow water analogue of SASI 
In view of the similarity of the equations, the formalism 
of advected-acoustic cycles studied by Foglizzo (2001, 
2002, 2009) can be translated into a linear coupling pro¬ 
cess between surface gravity waves and vorticity per¬ 
turbations, to produce a cycle between the hydraulic 
jump and the inner boundary. How would the absence 
of feedback from the advection of entropy perturbation 
affect the stability of this cycle in the shallow water 
analogue? The perturbative study and numerical sim¬ 
ulations of the shallow water equations performed by 
Foglizzo et al. (2012) confirmed the existence of an in¬ 
stability very similar to SASI, named SWASI as a shal¬ 
low water analogue of a shock instability. Despite the 
simplicity of the shallow water approximation, an excel¬ 
lent agreement was found between the oscillation period 
measured in the experiment and the one predicted by 
the shallow water equations. 


6.2 Experimental limitations 

6.2.1 Input parameters 

The variable parameters of the experiment are the flow 
rate 0 < Q < 4L/s, the thickness of the injection slit 
0.5 < H inn < 2mm, the height —5.6cm < z ns < 0 and 
diameter 8 < D ns < 16cm of the inner cylinder. The ra¬ 
dius of injection is set to 33cm. The minimum size of 
the injection slit is limited by its building accuracy. The 
diameter of the inner cylinder sets the maximum flow 
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rate above which the inner boundary should be mod¬ 
eled as a drowned boundary rather than a critical point. 
A central spinner, made of a light object in free rota¬ 
tion along a central vertical axis, is in passive contact 
with the surface of the water in the vicinity of the inner 
cylinder. Its rotating motion helps visualize the angu¬ 
lar momentum in the vicinity of the neutron star. In 
the new version described in Sect. |6.3.4[ the rotation of 
the full experiment is motorized and can be adjusted 
electronically. 


6.2.2 Hydraulic jumps vs shock waves 
The hydraulic jump in the experiment is a reverse ver¬ 
sion of the familiar circular hydraulic jump observed 
in kitchen sinks. The water flowing inward faster than 
the speed of waves (Fr > 1) is abruptly decelerated to 
slower velocities (Fr < 1) into a thicker layer of water. 
Mass flux is exactly conserved across this transition, 
and momentum flux is approximately conserved due to 
the minor drag on the fountain surface. Jump condi¬ 
tions similar to the Rankine-Hugoniot conditions are 
deduced from these two conservation laws. The conser¬ 
vation of the energy flux is not expected since the radial 
flow of water is described by only two variables, v and 
H , without any equivalent of the entropy variable in a 
gas. The same is true for isothermal shocks in astro¬ 
physics, where the excess energy is usually assumed to 
be radiated away. In the hydraulic jump with a mod¬ 
erate Froude number, excess energy is viscously dissi¬ 
pated within the thickness of the jump by one or several 
horizontal stationary rollers. Other processes such as 
wave emission, the formation of bubbles or splashes can 
participate in the evacuation of the excess energy at a 
higher Froude number (e.g. Chanson 2009). The 2D for¬ 
mulation of the shallow water equations does not resolve 
the internal structure of the hydraulic jump and simply 
assumes that energy is instantly dissipated across an 
idealized discontinuity. 

The description of hydraulic jump as a discontinuity is 
also unable to describe the limit of Froude numbers ap¬ 
proaching unity 1 < Fr < 2, where the radial extent of 
the hydraulic jump increases with a succession of oscil¬ 
lations referred to as an undular hydraulic jump. Ac¬ 
cording to the classification of Chow (1973), the range 
of Froude numbers most favorable to the simplest shal¬ 
low water description of the experiment is 4 < Fr < 9. 


6.2.3 Viscosity effects 

The strongest effect of water viscosity takes place ahead 
of the hydraulic jump, in the shallowest region where 
the vertical shear is responsible for a significant vis¬ 
cous drag. In the experiment presented in Foglizzo et al. 
(2012), this viscous drag has been measured and mod- 
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eled as a laminar drag in Eq. (13): 

„2 


dv 

— + (V xw)xu + V 


v 

y +c 


$ 


= a^w^p(19) 


where v w ~ 0.01cm 2 s _1 is the kinematic viscosity of wa¬ 
ter and a ~ 3 is a dimensionless number measured in 
the experiment. 

Fortunately this viscous drag is very small in the deeper 
region after the hydraulic jump and does not signifi¬ 
cantly affect the dynamics of the interaction between 
surface waves and vorticity perturbations. In order to 
compensate for the viscous drag ahead of the hydraulic 
jump, the water can be injected from the outer bound¬ 
ary with a higher velocity than the local free fall ve¬ 
locity, so that the velocity immediately before the jump 
Rj p ~ 20cm coincides with the local inviscid free fall ve¬ 
locity vg(R ip ) = (2 gH^(Rjp)) 1 / 2 . 

Increasing the injection flow rate increases the devel¬ 
opment of turbulence in the fountain, with an ex¬ 
pected transition for a Reynolds number of the order 
of Re ~ 2000. The viscosity of water seems too weak to 
affect the dynamics of SASI in the laminar regime, but 
a turbulent viscosity could significantly affect the effi¬ 
ciency of the cycle between surface waves and vorticity 
perturbations. 


6.2.4 Uncertain vertical profile of velocity 

The inviscid formulation of shallow water equations cor¬ 
responds to an idealized situation where the velocity v is 
uniform in the vertical direction. The non-uniform verti¬ 
cal velocity profile viscously induced by the no-slip con¬ 
dition introduces a modification of the quadratic terms 
such as (V x v) x v and Vv 2 /2 in the shallow water 
equations, since the vertically averaged value of v 2 is 
no longer equal to the square of the vertical average 
of v. Correction coefficients known as the Coriolis and 
Boussinesq coefficients could account for this effect if 
the vertical profile of velocity were known, for exam¬ 
ple in the idealized cases of a laminar flow or in fully 
developed turbulence. 

6.2.5 A 2D view of the equatorial plane of the 
stellar core 

The shallow water experiment is intrinsically limited to 
2D since the radial attraction of the central neutron 
star is mimicked from a projection of the vertical grav¬ 
ity of the laboratory onto the inclined surface of the 
fountain. The shallow water equations could, however, 
be solved in 3D if necessary, for example to assess the 
importance of buoyancy by comparison with 3D simu¬ 
lations of SASI. 


6.2.6 From shallow water equations to the 
experiment and back 

From the point of view of understanding the essence of 
the physics of SASI using analogies, the idealized set of 


inviscid shallow water equations is the simplest formu¬ 
lation with the most direct astrophysical connections. 
These equations can be studied analytically and solved 
numerically without using the experimental fountain. 
The water experiment is more complex because of vis¬ 
cosity effects which are not trivially captured in models 
and have no direct astrophysical interpretation such as 
the development of a boundary layer, the viscous drag, 
Coriolis and Boussinesq coefficients, the radial exten¬ 
sion of the hydraulic jump and its connection to bubble 
entrainment, splashes or unsteady behavior unrelated 
to SASI. A careful interpretation of experimental results 
is needed to identify the results which are governed by 
the simplest shallow water equations. The experiment 
is complementary to simulations because its limitations 
are different from the numerical ones, such as the dif¬ 
ficulty of first order numerical convergence due to the 
presence of a shock (Sato et al. 2009). This difficulty is 
inherent to any shock capturing method. The simplic¬ 
ity of use of the experiment makes it very convenient 
to rapidly explore a large parameter space and identify 
new phenomena. From an astrophysical perspective, the 
immediate next step consists in plugging the experimen¬ 
tal parameters of this new phenomenon in a numerical 
simulation of the simplest shallow water equations with 
a viscous drag. If the phenomenon is confirmed numeri¬ 
cally, its sensitivity to the viscosity is tested by perform¬ 
ing new shallow water simulations in the inviscid limit. 
The next step could be a 2D equatorial simulation of 
adiabatic gas dynamics to incorporate buoyancy effects, 
and a series of simulations to successively incorporate 
a hard inner boundary with a cooling function, neu¬ 
trino driven convection with a heating function, a third 
dimension, an equation of state for nuclear matter, neu¬ 
trino transport and interactions, relativistic corrections, 
a progenitor profile for the initial conditions and an in¬ 
ner boundary condition which allows the contraction of 
the proto-neutron star. It is important to bear in mind 
the extreme simplicity of the shallow water setup com¬ 
pared to the most advanced simulations such as Hanke 
et al. (2013). It is equally important to bear in mind 
that even the simplest setup like the shallow water for¬ 
mulation contains interesting hydrodynamical puzzles 
and surprises. 


6.3 Building up our physical intuition about 
SASI 

6.3.1 SASI is not familiar yet 

Part of the apparent complexity in the physical prob¬ 
lem of core collapse comes from hydrodynamical insta¬ 
bilities such as neutrino driven convection and SASI, 
which break the spherical symmetry. Most physicists are 
familiar with the convective instability, even though its 
interaction with a shock wave and an advection flow is 
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Figure 5. The coupling between vorticity perturbations and sur¬ 
face gravity waves in shallow water can help understand the cou¬ 
pling between vorticity perturbations and acoustic waves in a gas. 
The vorticity perturbation is shown in green at three successive 
positions in the upper illustration, viewed from above. The change 
in water elevation zb 5H produced by the vortical motion over the 
gradient of depth is viewed horizontally in the lower illustration. 
It is a source of surface gravity waves. 


far from trivial. By contrast, SASI is much less familiar 
to astrophysicists. Its advective-acoustic mechanism re¬ 
lates to aeroacoustic instabilities in ramjets (Abouseif et 
al. 1984), rocket motors (Mettenleiter et al. 2000), and 
can be traced back to the whistling of a kettle (Chanaud 
& Powel 1965). In astrophysics the instability of Bondi- 
Hoyle-Lyttleton accretion also relates to this category of 
unstable advective-acoustic cycles (Foglizzo et al. 2005), 
but the conical geometry of the shock make it less in¬ 
tuitive than SASI and this problem has been much less 
studied than stellar core collapse. Characterizing the 
fundamental properties of SASI in their simplest for¬ 
mulation makes it possible to enlighten more complex 
simulations (e.g. Blondin & Mezzacappa 2007). 

6.3.2 The linear coupling between vorticity 
perturbations and waves 

The translation of gas dynamics into shallow water dy¬ 
namics is in many respects a fruitful exercise, which 
can produce physical pictures sometimes more intuitive 
with water than with a gas although they are abso¬ 
lutely equivalent in the end. Let us describe the cou¬ 
pling between vorticity waves and surface gravity waves 
in an inhomogeneous flow. A vortical perturbation can 
be advected with the flow at uniform velocity without 
ever affecting the pressure equilibrium. Let us imag¬ 
ine that this perturbation is slowly advected in a di¬ 
rection where the depth of the shallow layer of fluid 
increases, as illustrated in Fig. [5j Because we have an 


intuitive representation that the surface of slow water 
(i.e. Fr <C 1) remains horizontal even when the depth 
increases abruptly (e.g. dH^/dr 1), we are able to 
guess that the vortical motion must perturb this flat 
surface by transporting some water from the deep re¬ 
gions to shallower regions and conversely. We also have 
an intuitive understanding that the perturbed water 
surface will return to horizontal through the propaga¬ 
tion of surface waves. The same intuitive reasoning can 
be translated into an isentropic or isothermal gas decel¬ 
erated by some external potential. A similar approach 
was used by Foglizzo & Tagger (2000) to explain the 
feedback from the advection of an entropy perturba¬ 
tion across the adiabatic compression produced by an 
external potential. It should be noted that these two 
feedback processes are fundamentally different but since 
entropy and vorticity perturbations are advected with 
the same velocity, the feedback they produce adds up 
in a coherent advective-acoustic cycle studied analyti¬ 
cally by Foglizzo (2009) and numerically by Sato et al. 
(2009). If a magnetic field were present, vorticity and 
entropy perturbations would no longer propagate at the 
same velocity and the advective-acoustic cycle would 
be transformed into five distinct MHD cycles (Guilet & 
Foglizzo 2010). 

6.3.3 Open questions accessible to the SWA SI 
fountain 

The SWASI fountain is currently being used as a com¬ 
plementary tool to numerical simulations to address the 
following questions: 

(i) Comparing the growth rate of SWASI in the ex¬ 
periment to the growth rate deduced from the shallow 
water equations with various prescriptions for the vis¬ 
cous drag and the turbulent viscosity can teach us about 
the sensitivity of the SWASI mechanism to experimen¬ 
tal limitations. 

(ii) Understanding the interaction of SWASI with the 
turbulence produced by the vertical shear may shed 
light on the physics of SASI and its interaction with 
the turbulence induced by parasitic instabilities or by 
neutrino driven convection. 

(iii) Measuring the experimental conditions for the 
symmetry breaking between the clockwise and anti¬ 
clockwise modes of SWASI can help us understand the 
nature of this non-linear process. 

(iv) Understanding the saturation amplitude of SASI 
is essential to predict which type of progenitors could 
lead to the strongest SASI oscillations. Is the growth of 
parasitic instabilities proposed by Guilet et al. (2010) 
always the dominant saturation mechanism? Testing 
the saturation mechanism of SWASI in the shallow wa¬ 
ter setup may be simpler than SASI in a gas because 
of the action of a single parasitic instability (Kelvin- 
Helmholtz) associated with the fragmentation of vor¬ 
ticity waves. 
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Figure 6. In 2014 the SWASI experiment was presented by the 
SN2NS collaboration (here R. Kazeroni) in the Paris Science Mu¬ 
seum to explain supernova theory to the public. 


6.3.4 new experiment including the angular 
momentum of the stellar core 

A new experiment has been built at CEA Saclay with 
similar dimensions as the one published by Foglizzo et 
al. (2012). The main difference is the capability to mo¬ 
torize the rotation of the full experiment so that water 
can be injected with a non-zero angular momentum. 
This device should be able to experimentally address 
the effect of rotation on the growth of the prograde spi¬ 
ral mode (Yamasaki & Foglizzo 2008) and its saturation 
amplitude. Based on the adverse rotation of the central 
spinner observed by Foglizzo et al. (2012) in a flow with¬ 
out angular momentum, one can expect to experimen¬ 
tally demonstrate the competition between the angular 
momentum of the background flow and the adverse an¬ 
gular momentum in the region close to the neutron star 
(Blondin & Mezzacappa 2007). 

6.4 A tool for public outreach 

A simplified version of the SWASI experiment was ex¬ 
posed to the public of the Science Museum in Paris, 
from December 2013 to February 2014 (Fig. |6|. The 
fountain was accompanied by a series of explanatory 
posters, a video about cosmic nucleosynthesis and a 
bouncing ball experiment. Two to four presentations of 
40mn were proposed every day by the 12 researchers of 
the SN2NS collaboration for an average audience of 15 
people aged 7 and above. It is remarkable that some of 
the most recent discoveries in supernova physics could 
be transmitted to more than 2,000 visitors during this 
period. The interest of the public was clear enough for 
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the museum to decide to incorporate this supernova 
fountain into their permanent collection as from 2015. 
This supernova fountain was awarded the 2014 prize 
for scientific communication from the French Ministry 
of Research and Higher Education. 

6-4-1 From bouncing balls to fluid mechanics 
As an introduction to the fluid experiment, the pub¬ 
lic was invited to comment on a side experiment were 
two bouncing balls fall vertically, both guided by a wire 
attached between the ceiling and the ground. The di¬ 
ameter of the lower bouncing ball is bigger (~ 8cm) 
than the upper one (~ 5cm). Used with a single ball, 
this experiment first illustrates the question of elastic 
vs inelastic bounce. When the two balls are dropped 
simultaneously the spectacular transmission of energy 
and momentum from the lower ball to the upper one is 
always a surprise to the public witnessing the ejection of 
the upper ball to the ceiling. This experiment is a good 
introduction to the simplest and naive solution to the 
problem of core-collapse. The SWASI experiment takes 
the public two steps further, both in the direction of 
fluid mechanics and with the additional degree of free¬ 
dom allowed by transverse motions. The radial flow of 
water in the experiment is presented as a view of the 
dynamics in the equatorial plane of the stellar core. 

6-4-2 Familiarity with water and kitchen sink 
hydraulic jumps 

One key aspect of the public success of the SWASI foun¬ 
tain is the simplicity of the setup, and the familiarity of 
the public with water. The main action on the fountain 
was a change in flow rate from 0 to lL/s which showed 
the radial injection of water, and the formation of a cir¬ 
cular hydraulic jump at a radius of ~ 20cm. The com¬ 
parison with the circular hydraulic jump visible in every 
kitchen sink helps illustrate the analogy with shocks, 
and identify the subsonic and supersonic parts of the 
flow. The changes in composition are also introduced 
during this phase of stationary accretion: while the su¬ 
personic gas is made of iron nuclei, the post-shock gas 
is dissociated into protons, neutrons and electrons, and 
the neutron star is essentially made of neutrons with 
neutrinos diffusing out. It is important to introduce the 
neutrinos as a key actor of the supernova mechanism 
early enough. This is one of the difficulties of the pre¬ 
sentations since this particle is generally unknown to the 
public, and absent from the analogous experiment. The 
post-shock gas is described as dense enough to capture a 
fraction of the outgoing neutrino flux, but insufficiently 
to revive the shock when distributed spherically. Unlike 
the stable jump in a kitchen sink, the hydraulic jump 
in the supernova fountain moves in an irresistible man¬ 
ner. Its motion is as inevitable as the fall of a pen in 
vertical equilibrium on a table. Like the falling pen the 
initial symmetry of the equilibrium is broken in a ran- 
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dom direction. The breaking of the circular symmetry 
appears with the development of m = 1 sloshing oscilla¬ 
tions along a random direction with a 3s period. As the 
oscillation of the hydraulic jump reaches larger ampli¬ 
tude, a second symmetry breaking leads to the rotation 
of the spiral wave in a random direction, clockwise or 
anti-clockwise. The central spinner raises the curiosity 
of the public by rotating erratically in the direction op¬ 
posite to the hydraulic jump. 

6-4-3 Capturing astrophysical processes accessible 
on a human scale 

The notions of time-scale and length-scale are intro¬ 
duced with the fountain, applying the factor one million 
to scale down the neutron star radius, the shock radius, 
the size of the iron core and the size of a supergiant. 
The factor ^hundred in timescale is applied to speed 
up the observed dynamics of the hydraulic jump. As the 
neutron star is identified with the surface of the inner 
cylinder, the extreme properties of these compact stars 
are stressed: nuclear densities, extreme spin rates, mys¬ 
terious velocities up to lOOOkm/s. The scaling factors 
enable the public to capture representations of extreme 
astrophysical processes on a human scale. 

6-4-4 Public access to the physics of explosion 
threshold, kick and spin 

The connection between the dynamics of the fountain 
and stellar explosions is explained along three observa¬ 
tions: 

(i) the explosion threshold is understood by observ¬ 
ing the radial extent of the post-shock material which 
can intercept neutrinos in the direction of largest shock 
displacement. Neutrino absorption will be enhanced in 
this direction compared to the spherically symmetric 
case. 

(ii) the spectacular pulsar kicks up to lOOOkm/s can 
be reached if the explosion takes place in an asymmetric 
manner: this is a consequence of the conservation of 
linear momentum, which can easily be experienced by 
feeling the recoil when throwing a heavy object 

(iii) the significant spin up of the neutron star is illus¬ 
trated by the motion of the central spinner. The conser¬ 
vation of angular momentum is invoked to accept that 
the spinner and the hydraulic jump rotate in opposite 
directions, noting that the injected flow is purely radial. 

7 CONCLUSION 

The development of 3D numerical simulations over the 
last two years has allowed a more detailed compari¬ 
son between theoretical models and observations. The 
most recent results suggest that the general framework 
of neutrino driven asymmetric explosions is consistent 
with several observations when the neutrino luminos¬ 
ity is adjusted to trigger the explosion. Among them, 


the kick of neutron stars (from Sclieck et al. 2004 to 
Wongwathanarat et al. 2013) and the chemical mixing 
in the stellar envelope (from Kifonidis et al. 2006 to 
Wongwathanarat et al. 2014). The distribution of 44 Ti 
clumps in the supernova remnant Cassiopea A (Grefen- 
stette et al. 2014) seems to confirm the theoretical pre¬ 
diction by Wongwathanarat et al. (2013). On the down¬ 
side, we should remember that we do not yet understand 
how to produce a robust supernova explosion from first 
principles. The last five years have shown several signs 
of convergence between a growing number of research 
groups. Part of the difficulties may come from the cost 
of checking the numerical convergence of 3D simula¬ 
tions in a dynamical system involving stochasticity and 
shocks in addition to the sensitivity to different numer¬ 
ical approximations. 

The dependence of the explosion mechanism on the pre¬ 
collapse structure of the stellar core has opened up the 
parameter space of initial conditions and stressed the 
reliance on the implementations of non-spherical effects 
in stellar evolution codes such as rotation, magnetic 
fields and convective asymmetries. The additional ef¬ 
fect of binary interaction during the lifetime of a mas¬ 
sive star may open the parameter space even more. The 
choice of initial conditions for 3D numerical simulations 
must be guided by a better understanding of the domain 
of dominance of the different liydrodynamical instabil¬ 
ities at work. Similarly the only way to assess if the re¬ 
sults from a simulation are general or coincidental is the 
detailed understanding of the mechanisms ruling their 
evolution. The shallow water analogy was developed as 
a complementary tool to improve our intuition about 
the dynamics of SASI. Its potential for public outreach 
has already been tested on the public with encouraging 
results. 

In the near future, supernova experts are impatient to 
discover possible clues through the properties of the 
compact object hiding behind the asymmetric ejecta 
of SN1987A. Even more direct constraints could come 
from the detection of neutrinos and gravitational waves 
from the next galactic supernova. There is ample work 
for theoreticians before these much awaited events oc¬ 
cur. 
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